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ABSTRACT 

A  control  thermodynamic  model  for  a  two-stage  Pulse  Tube  Refrigerators  (PTRs)  is 
developed  based  on  exergy  flow  in  the  refrigerator.  The  model  includes  flow  conductance, 
heat  transfer  effectiveness,  and  conduction  heat  transfer  parameters  for  regenerators  in  both 
stages.  It  is  assumed  that  a  phase  shift  controller  exists  in  both  stages  to  control  the  phase 
shift  between  the  pressure  and  mass  flow  rates  that  can  control  and  optimize  the 
performance  of  the  refrigerator.  The  effects  of  the  allocation  of  the  values  of  flow 
conductance  and  ineffectiveness  parameters  in  the  regenerators,  the  mid-stage  temperature, 
and  the  phase  shift  in  each  stage  on  the  performance  of  the  refrigerator  are  investigated. 
Important  dimensionless  parameters  controlling  the  thermodynamic  performance  of  two- 
stage  PTRs,  including  a  model  to  quantify  the  magnitude  of  different  components  of 
irreversibility  in  the  regenerators,  is  developed  and  discussed. 


KEYWORDS:  Exergy  analysis,  Multi-stage,  Cryocoolers,  Pulse  tubes,  Irreversibility 
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INTRODUCTION 


Multi-stage  Pulse  Tube  Refrigerators  (MPTRs)  are  attractive  cooling  systems  for  a  variety 
of  applications  and  their  staging  configurations  provide  the  flexibility  to  conveniently 
change  the  various  stage  temperatures  for  different  system  applications  [1-4].  Recently, 
interest  in  the  design  and  development  of  high  frequency  multistage  PTRs  has  been 
initiated  and  low,  no-load  temperatures  have  been  obtained  [5].  Numerical  calculation  to 
evaluate  the  performance  of  multi-stage  PTRs  has  also  been  reported  [6]. 

In  general,  the  model  of  control  thermodynamics  (finite  time  thermodynamics)  can 
be  applied  for  analysis  and  optimization  of  multi-stage  refrigerators  [7,  8].  Different 
thermodynamic  models  to  analyze  and  optimize  multi-stage  refrigerators  exist  in  the 
literature  and  use  different  models  to  calculate  the  irreversibility  for  internal  processes  and 
external  heat  transfer  with  reservoirs  [10-14],  to  just  name  a  few.  Thermodynamic 
optimization  of  two-stage  PTRs  has  been  carried  out  assuming  that  the  only  irreversibility 
in  the  regenerator  is  due  to  conduction  heat  transfer  [15].  In  this  study,  we  assume  that  a 
general  phase-shifter  exists  to  control  the  phase  shift  between  the  mass  flow  rate  and 
pressure  in  the  pulse  tube  for  each  stage.  In  addition,  we  include  a  more  realistic  model  of 
the  regenerator  in  both  stages  to  include  the  effects  of  the  void,  pressure  drop,  heat  transfer 
ineffectiveness,  and  conduction  heat  transfer  in  the  regenerator.  At  least  a  first  order 
thermodynamic  model  for  the  regenerator  is  necessary  to  be  able  to  provide  a  more  realistic 
thermodynamic  model  to  account  for  its  contribution  to  the  irreversibility  of  two-stage 
PTRs.  In  this  study  an  exergy-based  analysis  for  two-stage  PTRs  is  performed  to  quantify 
the  losses  and  optimize  the  refrigerator  [16-18].  Important  dimensionless  parameters 
controlling  the  performance  of  the  refrigerators  are  obtained  and  discussed.  In  addition,  the 
method  is  used  to  analyze  and  compare  two  generic  configurations  of  two-stage  PTRs 
including  the  effect  of  cooling  load  at  the  mid-stage. 


MATHEMATICAL  FORMULATION 

FIGURE  1  shows  a  generic  configuration  of  two-stage  PTRs  assuming  conductive 
cooling  at  the  midstage  using  two  separate  compressors.  This  generic  model  is  convenient 
for  thermodynamic  analysis  and  optimization  and  shows  the  important  characteristics  of 
two-stage  PTRs.  Exergy  comes  into  the  system  from  the  compressors  and  is  destroyed  as  it 
moves  into  the  system.  The  product  exergy  is  delivered  to  the  cold  reservoirs  at  the  mid¬ 
stage  and  the  second  stage.  If  there  is  no  cooling  load  at  the  midstage,  the  auxiliary  cooler 
only  supports  the  thermal  losses  from  the  first  regenerator.  We  assume  that  the  pressure 
and  mass  flow  rate  at  any  location  at  the  inlet  and  exit  of  each  component  are  given  by 

M  j  =  rhjCosioot  +  0 ;)  (1) 

P  j  =  Pa  +  PjCos(cot  +  6  j)  (2) 

where  0 .  and  0 .  are  the  phase  angle  for  mass  flow  and  pressure  at  any  inlet  or  exit  of  any 

component  in  the  system,  respectively.  The  pressure  component  of  exergy  at  any  location 
can  be  written  as  [18] 

RT  r  P 

<Ei>P~lMj\n(f)dt  (3) 
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Qc3>  TC3>  Pc3>™-c3 


FIGURE  1.  Schematic  of  the  two-stage  PTR  and  the  parameters  used  in  the  model 

In  most  practical  PTR  applications,  where  the  amplitude  ratio  is  relatively  small  (less  than 
0.3),  Equation  (3)  can  be  simplified  and  the  pressure  component  of  exergy  can  be  written 
as 

<  Ej  >,  -  (1/2 )TR0  rhj  ^  Cos(<pj  -0,)  (4) 

Pa 

The  thermal  exergy  transfer  at  the  cold  side  of  the  nth  regenerator,  due  to  conduction  heat 
transfer  and  ineffectiveness  of  the  regenerator,  can  be  written  as 

<  >lh=  (Qreg  +  Qcondn  X1  “  ~)  (5) 

where  Qregn  can  be  estimated  using  the  ineffectiveness  of  the  nth  regenerator  defined  by 

[18], 


Qre 


Qre 


A H„  (Thn  mhn  -  Tcn  mcn ) 


(6) 
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where  the  subscripts  hn  and  cn  denote  the  hot  side  and  cold  side  of  the  nth  regenerator  and 
Awis  its  ineffectiveness.  The  denominator  of  Equation  (6)  represents  the  maximum  rate  of 
enthalpy  transfer  into  the  regenerator.  Conduction  heat  transfer  to  the  cold  heat  exchanger 
can  be  estimated  assuming  a  linear  temperature  profile  in  the  regenerator, 

Qcondn  =  (KA/ L)n(Thn  -Tcn)  (7) 

where  (KA/ L)n  is  the  effective  thermal  conductance  of  the  nth  regenerator  and  its  shell.  It 
should  be  pointed  out  that  the  thermal  exergy  transfer  at  the  cold  side  of  the  regenerator  is 
opposite  to  the  direction  of  energy  transfer  at  the  point.  Using  Equations  (5)  to  (7)  the  total 
exergy  at  the  cold  side  of  nth  regenerator  can  be  written  as 

<  Ecn  >=<  Ecn  >p  +  <  Ecn  >,h~(m)RT0  mcn  ^ Cos{<j>cn  -  6J 

r  a 

+  [(l/ji)\nCp(Thn  mhn  -  Tcn  mj  +  (1 KA/L)n(Thn  -  Tj]  (1  -  U 

1  cn 

Oscillating  flow  of  heat  and  mass  transfer  in  regenerators  is  very  complicated, 
thermodynamic  analysis  we  consider  a  simple  model  for  the  regenerator  to  find  the  bounds 
for  cooling  capacity  and  efficiency  of  two-stage  PTRs.  It  is  assumed  that  the  mass  flow 
rate  at  the  cold  side  of  the  regenerator  can  be  obtained  using  a  linear  relation  based  on  a 
given,  appropriately  average  flow  conductance. 


(8) 

For 


mn+lCos(cot  +  <pn+l )  =  Cn[pnCos(wt  +  dn)~  pn^Cos(wt  +  0„+1)]  (9) 


where  C„  is  the  average  flow  conductance  in  the  «,h  regenerator.  In  addition,  a  simple 
model  for  the  conservation  of  mass  in  the  regenerator  can  be  written  as, 


mnCos{wt  +  <pn)  =  mn+lCos(cot  +  0„+1)  -  Vn  co[pnSin{cot  +  dn)  +  pn+lSin(wt  +  0,I+1)]  (10) 


where  V  n  is  the  properly  averaged  effective  regenerator  void  volume  including  the 
influence  of  temperature  distribution  in  the  nth  regenerator.  Expanding  Equations  (9)  and 
(10)  and  applying  trigonometric  identities  gives  a  set  of  equations  that  uniquely  determine 
the  relations  between  the  phase  shifts  and  amplitudes  of  pressure  and  mass  flow  rates  at  the 
cold  and  hot  sides  of  the  regenerator. 

Assuming  a  phase-shift  control  mechanism  exists  to  control  the  phase  shift  between 
the  mass  flow  and  pressure  at  the  cold  side  of  pulse  tubes,  one  can  optimize  the  exergy 
delivered  to  the  reservoirs  at  the  mid-stage  and  the  cold  stage  of  the  two-stage  PTR.  The 
irreversibility  in  the  pulse  tube  can  be  presented  by  defining  its  exergetic  efficiency  as  the 
ratio  of  exergy  transfer  at  the  hot  side  of  the  pulse  tube  to  the  cold  side  of  the  pulse  tube. 
Referring  to  FIGURE  1  and  using  exergy  balance  for  the  cold  heat  exchanger  at  the  mid¬ 
stage  of  the  two-stage  PTR  gives, 


<ec,  >=a,(U-  i)+-L 

Ui  *7 pt\ 


Qci+(\ht)\Cp(T0mhl -Tclmel)  +  (KA/L\(T0 -Tj 
+  (1  l*)hCf  (I,  mh2  -  T,  mc2  )  +  (KA/L)2  (Ta  -  Tcl ) 


(11) 
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where  subscripts  1  and  2  correspond  to  regenerators  1  and  2,  respectively.  Qcl  and  Tcl  are 
the  cooling  load  and  temperature  of  the  reservoir  at  the  midstage,  respectively.  The  bracket 
in  the  above  equation  represents  the  enthalpy  flow  in  the  pulse  tube  of  the  auxiliary  cooler 
and  rjpn  represents  its  exergetic  efficiency.  Equation  (8)  applied  to  the  first  regenerator 
gives  the  left  hand  side  of  Equation  (11).  Assuming  the  temperatures  of  the  cold  sides  of 
the  first  and  second  regenerators  are  the  same  and  equal  to  the  temperature  of  the  reservoir 
of  intermediate  cooling  load,  a  relation  for  the  cooling  load  at  the  mid-stage  can  be  written 
as, 


(\/2)RT.mcl  p  Cos((f>cl  -&cl) 

Qd  = - j - 2 — j - (l/x)A2Cp(7;wM  ~Tci  K 2 ) 

“  -1  + 

Tcl  Vpn 

- (KA! L)2 (Ta -Tcl  -(l/jt)\Cp(T0mM -Tclmcl)-(KA/ L\{Ta -Tcl) 


(12) 


where  Tcl  is  the  reservoir  temperature  of  intermediate  cooling  load.  Similarly,  application 
of  equation  (8)  to  the  cold  side  of  the  third  regenerator  in  the  second  stage  of  the  cooler  and 
the  balance  of  exergy  at  the  cold  heat  exchanger  of  the  second  stage  gives  a  relation  for  the 
cooling  load  at  the  second  stage, 


Qc2- 


(l/2)RT0mc3^Cos(<Pc3-dc3) 


- (Tcl  mc2  - Tc3  mci)- (KA/ L)3(Tcl -Tci) 


-1  + 


VpT2 


(13) 

Expansion  of  Equations  (9)  and  (10)  for  the  three  regenerators  gives  12  equations  relating 
the  amplitude  of  mass  flow  rates,  amplitude  of  the  pressure,  and  the  phase  angles  at  the 
inlet  and  exit  of  each  regenerator.  For  a  given  phase  shift  at  the  cold  sides  of  pulse  tubes, 
assumed  to  be  controlled  by  the  phase  shifters,  and  the  amplitude  of  pressure  at  the  driver 
sides  as  the  input  variables,  the  twelve  equations  can  be  used  to  find  the  amplitudes  of 
pressure  and  flow  rates  and  the  phase  shifts  at  all  other  locations.  Equation  (4)  can  be  used 
to  find  the  power  input  to  the  compressors, 


Wcompi  =  m)RTa  mkl  ^ Cos(<phl  -9J 

a 

(14) 

Kmp 2  =  (1  /  2 )RT„  mh2  CoS(<Ph2  -  dh2 ) 

(15) 

Using  Equations  (12)  to  (15),  the  performance  of  the  two  stage  pulse  tube  used  in  this  study 
can  be  evaluated.  The  exergetic  efficiency  of  the  two  stage  pulse  tube  can  be  calculated 
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Pressure  amplitude  ratio  across  the  second  stage  regenerator 
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FIGURE  2.  Exergy  delivered  at  the  second  stage  and  exergetic  efficiency  as  a  function  of  the  pressure 
amplitude  ratio  across  the  second  stage  regenerator 


from  equations  (12)  to  (15)  and  is  defined  by  the  exergy  delivered  to  the  cold  reservoir 
divided  by  the  total  compressor  power, 

QAT0/Tcl-l)  +  Qc2(To/Tc3-l) 


Vex  =: 


RESULTS  AND  DISCUSSIONS 


w  +w 

compl  comp2 


(16) 


Equations  (12)  to  (16)  give  the  basic  equations  for  the  calculations  of  the 
performance  of  the  two-stage  PTR  used  in  this  study.  The  mass  flow  rates  and  phase  shifts 
used  in  these  equations  are  obtained  from  the  trigonometric  relations  given  in  equations 
(14)  and  (15).  These  relations  are  considerably  simplified  for  the  case  of  no  void  in  the 
regenerator.  In  this  case  the  amplitude  of  mass  flow  rates  at  the  cold  and  hot  sides  of  each 
regenerator  are  the  same.  The  assumption  of  no  void  in  the  regenerators  provides  a 
thermodynamic  bound  for  the  performance  of  the  two-stage  PTR. 

There  are  several  parameters  that  can  be  controlled  or  used  as  variables  to  perform 
parametric  studies  and  optimization  of  the  two-stage  PTR  used  in  this  study.  Important 
dimensionless  parameters  that  results  from  equations  (9)  and  (10) 
area,,  =  Vn<a/Cn,Pn  =  mj Cnphn,mrn  =mcJmhn,  and  Pr„  =pcnl phn.  FIGURE  2  shows 

the  exergy  delivered  to  the  second  stage  reservoir  defined  by  Qcl  ( To  /  Tc03  - 1)  and  the 

exergetic  efficiency  of  the  two-stage  PTR  as  a  function  of  the  pressure  amplitude  ratio 
across  the  regenerator  of  the  second  stage.  In  this  calculation  no  load  is  applied  to  the  mid¬ 
stage  ( Qcl  =  0) .  In  the  calculation,  it  is  assumed  that  the  mass  flow  and  pressure  are  in- 
phase  at  the  cold  sides  of  both  pulse  tubes  and  the  values  of  ineffectiveness  for  the 
regenerators  are  fixed.  Important  values  of  parameters  used  in  the  calculation  are  given  in 
the  figure.  FIGURE  3  shows  the  cooling  capacity  and  the  exergetic  efficiency  as  a  function 
of  the  cold  end  temperature  of  the  second  stage  for  fixed  values  of  pressure  ratio  across 
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FIGURE  3.  Cooling  capacity  at  the  second  stage  and  exergetic  efficiency  as  a  function  of  the  cold  end 
temperature  of  the  second  stage. 


regenerator  3  and  the  mid-stage  temperature.  The  results  are  given  for  three  values  of 
second  stage  regenerator  ineffectiveness.  The  cooling  capacity  is  approximately  a  linear 
function  of  the  cold  reservoir  temperature  as  typically  occurs  in  practice.  The  rate  of 
increase  in  the  exergetic  efficiency  decreases  as  the  second  stage  temperature  approaches 
the  mid-stage  temperature.  FIGURE  4  shows  the  cooling  capacity  and  efficiency  diagram 
corresponding  to  FIGURE  2  when  the  pressure  ratio  across  regenerator  3  is  used  as  a 
parameter.  In  this  calculation  the  ineffectiveness  of  regenerator  3  is  used  as  a  parameter  to 
show  its  effect  on  the  performance  of  the  two-stage  PTR.  Other  parameters  are  given  in 
the  figure.  The  loop- shaped  curves  in  the  figure  indicate  a  compromise  between  the 
cooling  capacity  and  efficiency  of  the  two-stage  PTR  for  the  case  studied. 


cooling  capacity  at  the  second  stage  (W) 

FIGURE  4.  Cooling  capacity  and  efficiency  diagram  for  different  values  of  ineffectiveness  of  the  second 
stage  regenerator. 
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CONCLUSIONS 


A  control  thermodynamic  model  for  a  two-stage  PTR  is  formulated  using  the 
exergy  analysis  of  the  refrigerator.  A  simple  model  for  the  regenerators  was  developed 
based  on  given  flow  conductance.  An  important  dimensionless  parameter  for  parametric 
study  and  optimization  of  the  refrigerator  was  developed.  A  cooling  capacity  and 
efficiency  diagram  for  the  refrigerator  is  presented  and  the  effect  of  the  regenerator 
ineffectiveness  and  mid-stage  temperature  on  cooling  capacity  and  efficiency  is  evaluated. 
It  is  shown  that  depending  on  the  constraints  imposed  on  the  system,  a  compromise 
between  cooling  capacity  and  efficiency  is  possible.  Application  of  the  method  to  a  more 
conventional  two-stage  PTR  is  under  investigation. 
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